An explanation for the glass-like anomaly in the 
low-temperature specific heat of incommensurate phases 
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An explanation for the glass-like anomaly observed in the low-temperature specific heat of incom- 
mensurate phases is proposed. The key point of this explanation is the proper account for the phason 
damping when computing the thermodynamic magnitudes. The low-temperature specific heat of 
the incommensurate phases is discussed within three possible scenarios for the phason dynamics: no 
phason gap, static phason gap and a phason gap of dynamical origin. Existing NMR and inelastic 
scattering data indicate that these scenarios are possible in biphenyl, blue bronze K0.30M0O3 and 
BCPS respectively. Estimates of the corresponding low-temperature specific heat are in reasonable 
agreement with the experiments. 



Deviations from the Debye law observed in the low- 
temperature specific heat of incommensurate (IC) phases 
have attracted some attention 0, 0, ■ It is observed, 
in particular, a quasi- linear- in-T contribution to this spe- 
cific heat (C) at the lowest temperatures (< 1 K) and, in 
some compounds 0, a maximum in C(T)/T 3 at 
temperatures ~ 10 K. It is expected that these peculiar- 
ities are connected to the low-energy excitations specific 
to IC phases, i.e. the phasons. Theoretical discussions 
try to relate these phasons and the mentioned maximum 
111 C(T)/T 3 0,1. As a result, this maximum is usually 
ascribed to existence of a gap in the corresponding pha- 
son spectrum. The origin of the quasi-lincar-in-T contri- 
bution, in contrast, has not been elucidated. Neverthe- 
less, due to its similitude to that observed in glasses, this 
contribution inspire some speculations 0]. 

In this paper we show that these peculiarities in the 
low-temperature specific heat of IC phases have a natu- 
ral explanation. This explanation is based on a specific 
feature of phasons that is overlooked in the above cited 
works: the phason damping. In fact it follows from gen- 
eral considerations that, even in perfect crystals, phasons 
will be overdamped for small enough wavevectors 0]. 
The phason damping in perfect crystals has a strongly 
dependence on temperature and vanishes for T = (fj- 
Such a temperature-dependent damping, however, is not 
consistent with the NMR 0, 0] and the inelastic neutron 
scattering data 0, reported for the IC phases in 
which we arc interested: the damping inferred from these 
experiments is twelve orders of magnitude stronger than 
the calculated for perfect crystals (see below). The pha- 
son damping in real crystals may then be due to defects. 
To the best of our knowledge, there is no theoretical cal- 
culations of such a defect-induced phason damping. Any- 
way, to our purposes, this damping can be estimated from 
the above mentioned experiments. 

It is worth mentioning that, in some sense, phasons 



in IC phases are analogous to vortons in superconduct- 
ing vortex lattices. They both represent acoustical vi- 
brations of the corresponding IC structure that, unlike 
to the ordinary acustical vibrations, are overdamped for 
small enough wavevectors. In fact, it has already been 
found that vortons yield a linear-in-T contribution to the 
low-temperature specific heat of the corresponding vortex 
lattice 0, E| . The key point of the theory we present 
below coincides with that of the mentioned works on vor- 
tex lattices: the proper account for the specific thermo- 
dynamic features of the damped oscillators (phasons in 
our case). Mention also that these specific features have 
been found relevant even for "ordinary" crystals: optical 
vibrations in real crystals remain damped due to defects 
down to zero temperature, so they may give an impor- 
tant contribution to the corresponding low-temperature 
specific heat 01- I n this paper we shall concentrate on 
IC phases mainly because i) the smallness of the phason 
frequencies increases the importance of the phason damp- 
ing and, consequently, of the corresponding contribution 
to the specific heat and ii) the existence of experimental 
data characterizing the phason dynamics makes it possi- 
ble to estimate in order of magnitude. 

Let us first consider, for the sake of illustration, the 
simple case in which the phason spectrum shows an static 
phason gap (due to "phase pinning," see e.g. Ref. [l5j|). 
If the corresponding wavevector-dispersion is negligible, 
the phason branch becomes analogous to an optical one 
within the Einstein approximation. The corresponding 
contribution to the specific heat can then be computed 
from the specific heat of an harmonic oscillator C osc . But 
to describe correctly the low-temperature regime in the 
IC phases, it is crucial to take into account the phason 
damping. This can be done by using the formalism of 
Refs. LL2l Hal. As a result one obtains 
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where v = 2nk B T/h, ip'(x) = d 2 [lnT(x)]/dx 2 is the 
trigamma function, and the A's are the roots of the equa- 
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tion 



A 3 + u_dA 2 + ( w o + I^d)^ + WflW = 0, 



(2) 



with ujq the natural frequency of the oscillator and 7 the 
damping coefficient, i.e. the viscosity coefficient divided 
by the mass of the oscillator. Here lu^ 1 is a memory time 
associated with a Drude regularization of this viscosity 
damping (see, e.g., Ref. and the Appendix). In 

further considerations the limit lud —> 00 is taken. In 
this limit, the regularized viscosity damping reduces to 
the "ordinary" one. 

The dependence on temperature of the specific heat 
Eq. 0} of a damped harmonic oscillator can be seen in 
Fig. ^f 01 ' different values of the damping coefficient 7. 
In the undamped case (7 = 0), it is well known that 
this specific heat is exponentially small at low enough 
temperatures (fc^T <C fuvo). For any finite damping, 
however, the low-temperature asymptotic for the specific 
heat Eq. is linear-in-T: 



nk B 7 k B T 

T^Q 3 U>q hjjQ 



(3) 



In overdamped case (7 3> u)q) this expression is valid for 
fcgT <C fiwg/7 while in underdamped one (7 <C 1^0), it 
does for k B T <C hix>Q. 

It is worth noticing the qualitative similitude, already 
at this level of consideration, between this specific heat 
and that reported in IC phases 0, 0, 0, 13 I n particular 
in what concerns to the linear-in-T dependence at the 
lowest temperatures, which is connected to the damping 
as we have seen. In the following we shall make more 
concrete considerations. 

Let us then consider the Landau thermodynamic po- 
tential (see e.g. Ref. [tillO^ 

4> = ^ + ^ivl + vi) + ^if+vin |[(Vm) 2 + (v^) 2 ], 

(4) 

where rji and 772 are the real and imaginary part of the 
complex order parameter 77 = (171 , ^72) ■ In the IC phase 
(a < 0), the equilibrium values can be taken such that 



(eq) 



= (-a/b) 1 / 2 and rft 



(eq) 



0. Within the scheme 



of no phason gap, which seems to be valid for biphenyl 
(see below), for small deviations of the order-parameter 
components from their corresponding equilibrium values: 
7]{ = f?,- 00 ^ one has the following equations of motion 
(see e.g. Ref. @): 
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Hence 771 is associated with longitudinal (amplitude) fluc- 
tuations while r]2 does with transverse (phase) ones. 

In Eq. I|5b(l it is implicit that, as stated above, for 
small enough wave-vectors phasons are overdamped due 
to the viscosity term 7^7)2- At low enough temperatures 
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FIG. 1: Specific heat Eq. Q divided by T 3 of an harmonic 
oscillator for different dampings. Here units have been chosen 
such that fee = ft = 1, the natural frequency of the oscillator 
as ljq — 1 and the Drude frequency (see, e.g., Ref. [TH and 
the Appendix below) lud = 10 4 . A different scale in C osc /T 3 
is used in the plot for 7 = 10. 



the phason contribution to the specific heat will be lincar- 
in-T. From Eq. this contribution can be estimated 
as 
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for ksT/H <C mm(vk m , v 2 ^/^); where v 2 = c/m, 7 = 
7^/771 and k m represents the radius of the Brillouin zone. 
(Here it has been assumed that the damping is the same 
for whole phason branch.) This linear-in-T contribution 
to the low-temperature specific heat of IC phases due to 
phasons will prevails over the Debye one (due to acoustic 
phonons) at temperatures T < [j/(vk m )]9, where 9 = 
(H/kB)(V 3 /vY^km with V the velocity of sound. 

Let us mention that, in accordance with Eq. l|5af) . 
amplitude fluctuations of the order parameter (ampli- 
tudons) are also damped. In consequence, they also give 
a linear-in-T contribution to the specific heat at the low- 
est temperatures. However, this contribution is always 
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smaller than phason one because, for a given wavevec- 
tor, the amplitudon frequency is always greater than the 
phason one, the corresponding damping coefficient being 
similar in both cases [see Eqs. ||SJ| and ©]. Indeed the 
amplitudon contribution to the low-temperature specific 
heat of the IC phases of interest is not essential to de- 
scribe the corresponding experimental data, so we shall 
omit this contribution if further considerations. 

Let us now consider the scheme of a dynamic phason 
gap 0. In some compounds such as BCPS the coeffi- 
cients of the equations of motion for the order-parameter 
may have some frequency dispersion. A relatively simple 
dispersion is described in the phenomenological model 
developed in Ref. to reproduce the central peak in 
the high-temperature phase of strontium titanate. Such 
a dispersion can be understood as a result of the linear 
coupling between the order-parameter and a relaxational 
variable £ = (^1,^2) fiol fig. The phason dynamics is 
then governed by the following equations of motion: 

mi)' 2 + d 2 r]' 2 - cV 2 if 2 + d£' 2 = 0, (7a) 

7^2+^2+^2=0. (7b) 

Let us emphasize that the most important in these equa- 
tions is not the extra variable £2 by itself, but the 
frequency-dispersion that this variable yields in the pha- 
son dynamics. Strictly speaking, to reproduce the spe- 
cific features of e.g. BCPS, i.e. gapped phasons and a 
central peak, it is only this dispersion what would be 
necessary (see below). However, to calculate the corre- 
sponding phason contribution to the specific heat, it is 
convenient to bear in mind that such a dispersion is ob- 
tainable from Eqs. Q (see Appendix). Anyway, the 
resulting phason dynamic is different at different time 
scales: while the slow oscillations are overdamped (cen- 
tral peak), the rapid ones are not (gapped phasons) [Tgj . 
The phason contribution to the low-temperature specific 
heat then naturally split into two terms: C = C gp h + C cp 
(see Appendix for further details). The former: 
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is the contribution that can be ascribed to the (un- 
damped) gapped phasons 0]. Here wf (fc) = S 2 + v 2 k 2 , 
S 2 = d 2 /m, xq(T) = hS/(k B T) and, assuming that 
S <C vk m , O ~ Hvk m /k B . The latter contribution: 
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(9) 



is due to damped phasons (central peak). Here 
u; cp (k) = v 2 k 2 /(S 2 7i /m) and tf(x) = £[^hxT(x)] is 



the trigamma function. The low-T asymptotic of this 
latter contribution is given by Eq. 10 by replacing 
Jn 7€<5 2 - 

Prior to estimate the phason contribution to the low- 
temperature specific heat in particular cases, let us see 
why this scheme of a dynamic phason gap is plausible 
in BCPS. (In short, this was already mentioned in Ref. 
[nj.) From Eqs. Q one can find the susceptibility asso- 
ciated with the phason: 



X P h(k,uj)=m 1 [Z>l(k,uj) - to 2 - i~f{u>)w\ 



where uj 2 (k,ui) = uj 2 (k) 



(10) 



According to this susceptibility, the spectrum of the in- 
clastically scattered neutrons in the IC phase will show 
i) a central peak contribution and ii) a "gapped phason" 
contribution peaked at u> ~ 5 if <5 • Although the 

central peak is not resolved by neutron scattering [Tolllll| . 
the analysis of the corresponding intensities show that 
the situation in BCPS is compatible with this scenario 
[Toj . In what concerns to NMR experiments, one has to 
bear in mind that NMR is a low- frequency probe. The 
corresponding Larmor frequency is then expected to be 
such that ujl <C 7^ (it is actually the case in BCPS, see 
below). In consequence, the phason susceptibility probed 
in these experiments reduces to 



X P h(k,uj L ) ps m 1 (v 2 k 2 - lu 2 l - i^S 2 LU L ) 



(11) 



and, therefore, no gap will be revealed by NMR. (It can 
be said that only the central peak part of the spectrum of 
the inelastically scattered neutrons is probed by NMR.) 
It is worth mentioning that formula (|llfl was used in fact 

in Ref. @- This further gives a uo 1 ^ 2 dependence of the 
spin-lattice relaxation time, first derived in Ref. |2Cj . 

We are now in a position to estimate the phason con- 
tribution to the low-temperature specific heat in the IC 
phases of biphenyl and BCPS, and in the charge-density- 
wave system blue bronze K0.30M0O3: 

i) Biphenyl. In accordance with NMR and inelastic 
neutron scattering data , the scheme of no phason gap 
seems to be valid here (the corresponding gap should be 
< 10 MHz at least). From neutron scattering, the phason 
damping can be inferred such that 7 ~ 50 GHz. Notice 
that such a damping, obtained at T = 3 K, is twelve 
orders of magnitude stronger than the one expected for 
a perfect crystal . So defects must be the cause of the 
real damping, although inducing no (observable) phason 
gap. Both phason and sound velocities can be estimated 
as - 10 5 cm s" 1 H H, E3. Therefore, the linear-in-T 
contribution due to phasons will prevail over the Dcbyc 
one at temperatures < 1 K. In addition, according to 
Eq. ©, it is obtained C/T ~ 10 10 2 erg cm- 3 K" 2 
for these low temperatures. The experimental value ~ 
0.5 erg cm -3 K -2 is smaller in this case, what shall be 
commented below. 

ii) BCPS (dynamic phason gap). From NMR and in- 
elastic neutron scattering experiments H,EJEI> the cor- 
responding damping can be inferred such that 7^i5 2 ~ 
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10 GHz. Notice that no temperature dependence is ob- 
served in this damping down to T = 19 K, so it seems 
very probable that it remains unaltered down to very low 
temperatures (as in biphenyl) . According to Refs. [3.l2ll|. 
both phason and sound velocities are also ~ 10 cms -1 . 
Then, the linear-in-T contribution due to damped pha- 
sons is expected to be the leading one at temperatures < 
1 K. It would be such that C/T ~ 10-MO 2 erg cm" 3 K~ 2 , 
what agrees in order of magnitude with the observed ex- 
perimentally 0- In this case (dynamic phason gap), in 
addition to this linear-in-T contribution there is a contri- 
bution due to undamped gapped phasons [Eq. JS|)]. The 
phason gap can be estimated as 5 ~ 100 GHz [l(J, so the 
maximum in C/T 3 observed at ~ 1 2K can then be 
explained as a result of this latter contribution [l^ ■ This 
has already noticed in Ref. Q. 

hi) Blue bronze Kq^MoO^,. Here, as well as in other 
charge-density- wave systems (see, e.g., Ref. |22| and the 
references therein), the scheme of a static phason gap 
seems to be the most appropriate one. In accordance with 
inelastic scattering data [23j this gap is ~ 200 GHz, and 
the phason damping can be estimated as 7 ~ 800 GHz. 
Both phason |2J] and sound 0, 3 velocities can be in- 
ferred as ~ 3 x 10 5 cm s^ 1 . Consequently, a linear-in-T 
contribution with C/T ~ 10 -j- 10 2 ergcm~ 3 K~ 2 is ex- 
pected due to the damping of the phasons (prevailing 
over the Debye one at temperatures < 1 K) and, in addi- 
tion, a maximum in C/T 3 at T ~ 10 K due to the phason 
gap. This qualitatively describes the experimental obser- 
vations Ply)]. Mention that in Ref. j25j it is questioned 
the explanation that the maximum in C/T 3 can be as- 
cribed to gapped phasons: it is shown that a similar max- 
imum arises due to the acoustic anisotropy of this sys- 
tem. In Ref. [2(| it is pointed out that, however, acoustic 
anisotropy cannot explain why this maximum is so sen- 
sitive to the crystal quality It is also worth mentioning 
that, given the relatively large anisotropy in the phason 
velocity [24|, phasons can behave as overdamped oscilla- 
tions along certain directions. Bearing in mind that an 
overdamped phason branch in the one-dimensional case 
give rise to a contribution ~ T 1 / 2 in the low-temperature 
specific heat , this could explain the deviation of the 
linear-in-T behavior observed experimentally in Ref. Q 
at the lowest temperatures. 

It is worth mentioning that, at present, the above con- 
tributions to the low-temperature specific heat due to 
the damping of the phason can only be estimated in or- 
der of magnitude. It is true that the corresponding pha- 
son dynamics has been extensively studied in a number 
of papers. However, the quantitative comparison of our 
results with the experimental data would require, in par- 
ticular, the knowledge of the phason damping for the 
whole phason branch. Unfortunately, these latter data 
are not reported in the literature and consequently, when 
carrying out the comparison of our results with experi- 
mental data, we are forced to make some assumptions. 
For instance that the phason damping is the same for the 
whole phason branch. These type of assumptions, at the 



present unavoidable, could be the reason for the over- 
estimation of the linear-in-T contribution due the pha- 
son damping when comparing with experimental obser- 
vations. More precise estimates require more complete 
experimental data (see Ref. f° r a detailed discus- 

sion). Nevertheless, the agreement in order of magnitude 
we have obtained indicates that the deviations from the 
Debye law observed in the very-low-tcmpcrature specific 
heat of IC phases are, the most probably, in regard to the 
phason damping. Mention also that the resulting linear- 
in-T behavior is expected as a low-temperature asymp- 
totic limit. It is well possible that in the above mentioned 
experiments this limit has not been fully achieved but 
rather arc in an intermetiatc (crossover) region, which 
could be the reason of the different power laws observed 
there. 

In summary, we have shown that the glass-like anomaly 
observed in the low-temperature specific heat of incom- 
mensurate phases can be explained as a result of the 
phason damping. Three possible scenarios for the pha- 
son dynamic, reproducing the corresponding NMR and 
inelastic scattering data for biphenyl, BCPS and blue 
bronze K0.30M0O3, are discussed. Estimates of the corre- 
sponding low-temperature specific heat are in reasonable 
agreement with the observed experimentally 

We thank K. Biljakovic for useful discussions. 



APPENDIX A: PHASON CONTRIBUTION TO 
THE SPECIFIC HEAT 

The knowledge of the partition function Z of a system 
permits the computation of the thermodynamic magni- 
tudes. We are interested, in particular, in the specific 
heat: C = ^[d 2 (T In Z)/dT 2 ], where (3 = l/(k B T). 

Formally, Z can be written as the path integral [ihl 

Z = l Dfa(r, r)] exp (-£[&]/») , (Al) 

where 0j(r,r) (r € V, r £ [0, H/3]) arc the fields char- 
acterizing the system, S[<pi] is the Euclidean (imaginary- 
time) action of the system, and the path sum is over the 
fields satisfying the appropriate boundary conditions. In 
what concerns to the /i/3-space of the (imaginary) time 
variable r, these boundary conditions are periodic ones 
for fields associated with bosons (to our purposes, it is 
sufficient to consider here bosonic systems). 

Among the fields characterizing the system, it can be 
identified the one associated with the degrees of freedom 
of interest, say 77(1*, r). Formal integration over all other 
degrees of freedom, i.e., over the set of fields {c/5,(r,r)}, 
gives 

Z=l 2%(r, t)\ exp (-S eS [r)]/h) , (A2) 



where the effective Euclidean action S^fffa] is such that 

cxp (-S eS [n]/hj = <j> P[^(r,r)]exp (S[t], tpi]/h) . 

(A3) 

This action can be divided into S^sM = S^ff + Sefffa]- 
Thus the partition function can be expressed as Z = 
Z Z red , where Z Q = cxp ( - /ft) and 

Z K d = I V[rj(r, t)} exp (-S eS [r,}/K) . (A4) 

The reduced partition function Z lcd results from the con- 
tribution due to the subsystem characterized by ?y(r, r). 
In contrast to the action in Eq. (|A1|) , the effective action 
5 e ff [rj\ in Eq. i|A4[l may be, e.g., a temperature-dependent 
functional. It is because S c s[ri] arises from a partial inte- 
gration. This procedure is well known, e.g., in the theory 
of phase transitions (see, e.g., Ref. [3lT | for an exposure in 
the classical case). It is worth mentioning that Sos^] has 
not to be hcrmitian necessarily. Relaxational processes 
can then be described within this formalism. We are in- 
terested, in particular, in the contribution to the specific 
heat due to phasons (damped oscillators). Let us illus- 
trate the method that further we shall use to calculate 
such a contribution considering first the case of "normal" 
phonons (undamped oscillators). 

1. Noninteracting undamped phonons: 
Path-integral calculation of the partition function 

In the simplest case, the effective Euclidean action for 
phonons can be taken as 

S eS [v) = / *jf (j + \^ + frr,f) > 

(A5) 

where ry(r, r) is a (bosonic) displacement field. The first 
term of this action is associated with the kinetic energy 



while the last two terms do with the potential one. The 
case a = corresponds to acoustic phonons while a =/= 
does with optic ones. 

The first variation of the action gives the corresponding 
equation of motion. According to Eq. (| AB|> it reads 

mAr) = to (-8$ v 2 V 2 ) r) = (A6) 

in the ft/3-space, where uj 2 = a/m, v 2 = c/m, and A is 
the so-called fluctuation operator (the reason of such a 
designation can be seen below). 

By putting i] = r/o + rf in Eq. (|A5|) , with ?yo satisfying 
the corresponding equation of motion together with pe- 
riodic boundary conditions, and integrating by parts, we 
find that 

to f r hf3 

S eS [rj\ - S eS [r)o) + YJ dr J dT11 ' Av ' (A7) 

for periodic fluctuations ?y'. Consequently, the reduced 
partition function (|A4|) can be expressed as Z rcd = 
Z °ed Z ^ whcrc Z °ed = exp(-S eff M/ft) and Zf cd is as- 
sociated with fluctuations: 

Zf ed = ^©[V(r,r)]exp (~ J dr J™ drrfhA 

(A8) 

(Gaussian fluctuations in our case). 

The only solution of Eq. (|A6|) which is periodic with 
respect to the r variable reduces to rjo = 0. So that 
Z° cd = 1 in our case. Expressing rj' as the Fourier trans- 
form 

oo 

v'(r,r) = J2 E ^u^'e^, (A9) 

k=0 n= — oo 

where v n = 2im/ (ft/3) are the so-called bosonic Mat- 
subara frequencies |32| ; and writing down the functional 
measure as j^H 



V[ V '(t,t)] 



n 

k=0 



d vLo 



v /2 7 rft 2 /3/TO 



n 



Tr/(m/3v%) 



(A10) 
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we have 

Z rc d = J 



k=0 



dl l'k,o -q ff f dRe rfk^dkn Vk, v „ 




2 

k=0 n=0 



= n 2sinh[^ (fc)/2]' ^ Al1 ^ 

where A(fc,z/ n ) = + an< ^ = + u 2 fc 2 (note that the reduced partition function we obtain is exact 

because it is connected to a quadratic effective action) . In the last form, we have the usual expression for the partition 
function of a system of noninteracting undamped phonons (see, e.g., Ref. [23). 

From Eq. (|A6|) it can be seen that the eigenvalues of the operator A gives the inverse of the (linear) dynamic 
generalized susceptibility: x _1 (k, lj) = mA(k, u>). When acting on fi,/3-periodic functions, these eigenvalues can be 
obtained by the transformation 

A(k,w) — > A(k,ii/ n )=A(k,i/ n ). (A12) 

[In the case of phonons: A(k, uS) = u>o(k) — uj 2 — > A(k, iv n ) = A(k, v n ) = ui 2 (k) + v 2 .] As we have seen, the 
knowledge of these latter eigenvalues permits the computation of the reduced partition function. It is the case even 
for relaxational dynamics. In the following, we shall take advantage of this point when calculating the contribution 
to the specific heat associated with phasons. 

2. Phason contribution to the specific heat 

a. Scheme of no phason gap 

Let us now calculate the phason contribution to the specific heat within the scheme of no phason gap. The procedure 
we shall follow is based on the results concerning to the damped harmonic oscillator, which can be found, e.g., in 
Rcfs. ^|[2l|. Formulas we generate are analogous to that given in Ref. [T3 | . 

According to Eq. I|5bjl the phason susceptibility is 

X P h(k,uj) = 1 : — =-, (A13) 

m\U>2(K) — LO z — 

where Lo 2 (k) = v 2 k 2 , v 2 = c/m, and 7 = 7,,/m. From this susceptibility one finds that 

A ph (fc,w) = [mxphik,^)}' 1 

_ [lu + i\i(k)][u + i\ 2 (k)}[uj + i\ 3 (k)} 



where 



(A14) 



\i{k) + \ 2 {k) + \ 3 (k) =uj d , (Al5a) 
Ai(fc)A 2 (fc) + A 2 (fc)A 3 (fc) + Ai(ft)As(*) = u$(k) + >yu D , (A15b) 
Ai(fc)A 2 (fc)A 3 (fc) = uj D uj 2 (k), (A15c) 



and uj d is a memory time associated with a Drude regularization of the damping term in Eq. J5bj): 

7 v m(t)^i v UD [ dt'e-^-^mit'). (A16) 

J —00 

Below we shall take the limit u>d —> 00 in which this regularized damping term coincides with the original one. But, 
at this start point, such a regularization is necessary in order that the infini te-p roduct expression for the partition 
function of the system, see Eq. (|A17(I below, actually exists (see Ref. [l3llT6ll3^ ). 
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The phason contribution to the partition function of the system can be computed from A p h by going to the space 
of fi/3-periodic functions, i.e., by replacing lu —> iv„ with v n = 2irn/(Kf3) [Til l29j (as we have seen, A then gives the 
eigenvalues of the so-called fluctuation operator [lfij): 

1 °°2 \ / 1 °° 2 / \ \ 



z ? h - II ( r a / a . n k\ II rr irzr\ ) - II ( i^zm II 



k= x \npy/ A ph (k,o) A ph (fc, iu n ) J ii ynpu 2 (k) i = i K + Ai(fc)]K + A 2 (fc)]K + M*)] / 
- ii [ r^) J ' (A17) 

where f = v\ and r(x) is the gamma function. We then have that the corresponding contribution to the specific heat 



is 



C = k B T d2(T ^ 2 Zph) =k B J ^j^^WHVtA^)/^] - (^A)VKM - l}, (A18) 

where ip'(x) = d 2 [hiT(x)]/dx 2 is the trigamma function and summation over wavevectors has been replaced by 
integration. 

By inserting the asymptotic expression x 2 ip'(x) rs x + 1/2 + l/(6x) into the integrand in Eq. I|A18(1 and using 

x — >oo 

the relations (|A15() . we further get the asymptotic low-temperature specific heat due to phasons: 



1,3 L2 T 

C « ^ITWTv (A19) 
t^o 6?r H(v 2 k 2 Jj) 



b. Scheme of dynamic phason gap 

Let us finally calculate the phason contribution to the specific heat within the scheme of a dynamic phason gap. 
Instead of Eq. (|7b|) , let us consider the equation 

l&W + f dt' ls (t - t%(t') + &(*) + d m (t) = 0, (A20) 

where J^(t) = "i^lod exp(— ujut). Here we have introduced a inertial term and considered a Drude regularization for 
the damping one. Similar to the case of no phason gap, below we take the limit (i — > and u>d — ► oo in which Eq. 
(|A20|) reduces to Eq. ffb)l . 

From the corresponding susceptibilities we get 



A ph (fc,w) = m^[x m \ 2 (k,uj)x i2 \ 2 (k 1 uj) - X V2 %(k,^)] = [wf (fc) - lj 2 ](w 2 - lu 2 - iw 2 ) ~ S 

_ [uj + i\i{k)][uj + i\ 2 {k)}[m + iX 3 (k)][uj + iX 4 (k)][uj + iX 5 (k)] 

uj + icon 



2 2 



(A21) 



where w| = S 2 + v 2 k 2 , w 2 = and now 7 = j^/m. 

Within the scheme of the dynamic phason gap we are considering, the phason contribution to the partition function 
of the system must be understood as the contribution of the degrees of freedom represented by ?/ 2 an d £,2, i-c, the 
contribution of linearly coupled oscillators. Following Ref. 0,|2!| one can see that it can be computed from A p h: 



^=n — ^— ft 



=0 V(fy3) 2 ^/A ph (fc,0) «=1 A ph(M^ 



k= \ {H0) 2 wy/v2(k) - 6 2 ^ K + Ai(fe)][i/„ + A 2 (fc)]K + A 3 (fc)]K + A 4 (A;)][^ + A 5 (fc)] J 

( (to\2 I Ai(fc)A2(fc)A 3 (fc)A 4 (fc)A 5 (fc) N \ 1/2 

tt m { •*» ) r[A 1 (fc)/Hr[A 2 (fc)/Hr[A 3 (fc)/Hr[A 4 (fc)/Hr[A 5 (fc)/H , 

11 W fM I ' (A22) 
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We then have that the phason contribution to the specific heat is 



C = k B T 



a 2 (rinz ph ) 

g T 2 



d 3 k 
(2tt) 3 



J2[\i(k)M 2 i>'[\i(k)/v} - {u> D /v) 2 ^'{u> D /u) - 2 



i=l 

We are interested in the limiting case w 2 (fc) — 8 2 ,j 2 , u> 2 {k) <C w 2 <C oj 2 d - We then have 

Ai(fc) u> cp (k), A 2 (fe)-7W 2 -w cp (fc), A 3 w w D - ~fw 2 , A 4 , 5 (fc) ~ ±iw 2 (A;), 



(A23) 



(A24) 



where w C p(fc) = v 2 k 2 /(S 2 ^). A a result, in Eq. I|A23(1 the contribution due to the central peak (damped phasons) and 
the one due to the (undamped) gapped phasons naturally split up: C = C cp + C gp h- The low-T asymptotic of the 
former contribution coincides with Eq. I|A19() by replacing 7„ — > 7^<5 2 . The latter is given by 



Cgph — kg 
pa -k 

= ks 



<7 3 k r -1 

— {[Ar^/HVlAiW/H + [Mk)M 2 nHk)M 1} 
|^{i + 2M^] 2 Re n^(k)M} = k B J (™ 2 ^ 

kl fT\ 3 / >e/T J [x 2 - x 2 (T)] 1 / 2 x' i e x 



- 1 



2tt 2 ve 



dx- 



x (T) 



(e* - If 



(A25) 



where xq(T) = M/(fcsT) and O ~ Hvk m /kB (here it has been assumed that 8 vk m ). 
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— V-k,-v n - 

One could expect that such a measure merely be 



n 



n=l ^ J J 



dRe-n' kil , n dlmri' k ^ r 



The origin of the nontrivial factor appearing in Eq. HAIOII can be understood following, e.g., Ref. |3Cjl . The starting point 
should be, in fact, to consider the effective Euclidean action 



S eft [rj, 7r] = dr 



a/3 



in(d T n) 



instead of that given in Eq. 1A5I . Thus, in addttion to the field rj(r, r) associated with the displacement amplitudes, we 
take into account the one representing its conjugated momentum density tt(t,t). Since this action is quadratic in n(r, r), 
we can complete the square and further evaluate the path integral over 7r(r, r) in the reduced partition function: 



V[ V (r, T)]V[7v(r, r)\ exp -S ef t[rj, n]/H 



£>[77(r,T)]exp(-Seff[7?]A), 
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where SW^] coincides with 1A5|I and the new measure V[n(r, r)] is the old one T>[n(r, r)] multiplied by mentioned factor. 
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